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This report summarizes my recent work on the construction of finite- 
difference models of differential equations having zero truncation errors. 
The details of the calculations which follow will not be given since a paper 
is being written on these topics. 


I. Unidirectional 'lave Equation 

The unidirectional, linear wave equation in one-dimension is 

(1) u + u* = Xu 

X o 

X = constant. 

For the initial value problem:, where 

(2) u(x,o) = f(x) = given function, 
the exact, general solution is 

(3) u( x, t ) = e'^ f(x-t) 

If we define 

(4) u(x n ,t n ) = u" 

x = (Ax)m , t = (At)n 
rn n 


where (m,n) are integers, then it is easy to prove that the following finite- 
difference scheme is a zero truncation error model of ea. (l): 


, r s n+1 Xh n 
(5a) u m = e u, 


m-1 


Ax = At = h 


where 

(5b) u° = f ( (Ax) m ) 
m 
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It should be pointed out that the conventional use of either the 
Euler forward- or backward- difference schemes for either of the derivatives 
in eq. (1) does not lead to eq. (5a). In more detail, the ’’full” difference 
scheme associated with eq. (1) which has zero truncation error is easily 
gotten from eq. (5a) and written out in detail is 


I'M A^*' u n u n 

(6) t m + t m 


- Xu 


n 

m-1 


where 


(7a) A+ + ^ u" = u“ +1 - u‘ 


(7b) A v ; u 
x 


( +) n n 

'•'ll =u - 


n- i 


Note that the conventional 
eo. (1) would give 


mention or Unite amerence teermeu es 


( 8 ) 


A (t) r. 
A t u n 


A xW u ! 


At 


Ax 


or some similar result, which aces have a non-zero truncation error. Direct 
comparison cf eqs. (6) a^d (8) shows that previous knowledge of finite-difference 
models for differential equations may not be of great help in the construction 


of schemes with zero-truncation error. ’ 

For practical problems where eq. (1) is to be numerically integrated, 
either eq. (5a) or eq. (6) may be used since they are mathematical equivalent 
expressions. However, from a computational point of view, eq. (5a) should 


be used. 
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Recently, I have beer. able to obtain a zero truncation error 
model for the following nonlinear, uni-directional wave equation 

(9) u t + u v = X u'" (u = positive integer) 

The general cast is rather difficult to write down ; however, for n = 

(10a) Uj. u _ = \.r 

o X 

ana the difference scheme is 

n 

> * 



difference 


z, we have 


where 


